Spin torque switching of an in-plane magnetized system in a thermally activated 

region 
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The current dependence of the exponent of the spin torque switching rate of an in-plane mag- 
netized system was investigated by solving the Fokker-Planck equation with low temperature and 
small damping and current approximations. We derived the analytical expressions of the critical 
currents, I c and I*. At I c , the initial state parallel to the easy axis becomes unstable while at 
I*(— 1.27/c) the switching occurs without the thermal fluctuation. The current dependence of the 
exponent of the switching rate is well described by (1 — I/I*) b , where the value of the exponent b 
is approximately unity for I < I c while b rapidly increases up to ~ 2.2 with increasing current for 
I c </</*• The linear dependence for I < I c agrees with the other works, while the nonlinear 
dependence for J c < / < I£ was newly found by the present work. The nonlinear dependence is 
important for analysis of the experimental results, because most experiments are performed in the 
current region of I c < I < I*. 

PACS numbers: 75.78.-n, 85.75.-d, 75.60.Jk, 05.40.Jc 



I. INTRODUCTION 

Spin torque switching 1,2 of a magnetization in nanos- 
tructured ferromagnets has attracted much attention due 
to its potential application to spintronics devices, such as 
spin random access memory (Spin RAM)^. An accurate 
estimation of the thermal stability, Ao, is very important 
for these applications, where the thermal stability is de- 
fined as the ratio of the magnetic anisotropy energy to 
the temperature (fee? 1 )- For example, the retention time 
of the Spin RAM exponentially depends on the thermal 
stability^. 

The thermal stability is experimentally determined 
by measuring the spin torque switching in the ther- 
mally activated region^, and analyzing the switching 
probability^ with the formula P = 1 — exp[— dt'v(t')\, 
where v = /exp(— A) and / are the switching rate and 
the attempt frequency, respectively. The exponent of the 
switching rate, 



A 1 



(1) 



is the energy barrier height of the switching divided by 
the temperature. In this paper, we call A the switching 
barrier. The thermal stability is defined as Ao. Here, 
I, and I c are the current and the critical current of the 
spin torque switching at zero temperature. Equation ([1]) 
characterizes the switching in the thermally activated re- 
gion, is and applicable for |/| < \I C \, The exponent of 
the term 1 — I/I c in Eq. ([T]) is denoted as b. Equation 
([1]) was analytically derived for the uniaxially anisotropic 
system by solving the Fokker-Planck equation 8 - - — , and 
6 = 2. Recently, it has also been confirmed by numerical 
simulations^—. 



It is difficult to derive the analytical formula of the 
switching barrier for an in-plane magnetized system, 
which does not have axial symmetry due to the presence 
of two anisotropic axes (an easy axis along the in-plane 
and a hard axis normal to the plane). Previous analy- 
ses adopted b = 1, which has been obtained by solving 
the Landau-Lifshitz-Gilbert (LLG) equation 1 ^ as well as 
the Fokker-Planck equation^. The point of these works 
is that the effect of the spin torque on the switching 
barrier can be described by the effective temperature^ 
T/(l — I/Ic)- However, the effective temperature ap- 
proach is valid only for I <C I^^. The current depen- 
dence of the switching barrier in a relatively large cur- 
rent region (/ ~ I c ) remains unclear, while such large 
current is applied to a ferromagnet in experiments to 
quickly observe the switching. These facts motivated us 
to study the switching barrier in the large current re- 
gion. In RefJ^, we numerically solved the LLG equa- 
tion of the in-plane magnetized system, and found that 
the switching does not occur even if I > I c , although I c 
has been assumed to be the critical current of the spin 
torque switching at zero temperature. We also found 
that the exponent b is larger than 2 by analyzing the 
numerical results with a phenomenological model of the 
switchin g 17 : 18 . 

In this paper, we develop an analytical theory of 
the spin torque switching of an in-plane magnetized 
system based on the Fokker-Planck theory with WKB 



approximation: 



,19-23 



where the temperature is assumed 



to be low, and the magnitudes of the damping and the 
spin torque are assumed to be small. We derive the an- 
alytical expressions of the critical currents, I c and I*. 
The initial state parallel to the easy axis becomes un- 
stable at I — I c , as shown in Ref. 24 . On the other 



2 



hand, at I = 1.27/ c ), the switching occurs with- 

out the thermal fluctuation. We also find that the cur- 
rent dependence of the switching barrier is well described 
by (1 — I/I*) b , where the value of the exponent b is 
approximately unity for the current I < I c , while b 
rapidly increases up to ~ 2.2 with increasing current for 
I c < I < I*, showing a good agreement with Refsi 15 ' 16 . 
The nonlinear dependence for I c < / < I* newly found in 
this paper is important to evaluate the thermal stability 
because most experiments are performed in the current 
region of I c < I < I* . 

The paper is organized as follows. In Sec. [TT1 we 
summarize the Fokker-Planck equation of the magneti- 
zation switching. In Sec. IIII| we describe the details of 
the WKB approximation. Equations ((2T|) , ([19]) and (f20| . 
lead to the main results shown in the next sections. In 
Sec. IIV1 we apply the above formula to the uniaxially 
anisotropic system, and show that the present formula 
reproduces the results 6 = 2 derived in Refs^ ' 10 ' 11 . In 
Sec. [Vj the switching barrier of an in-plane magnetized 
system is calculated. In Sec. IVI1 we compare the current 
results with our previous works. Section IViTl is devoted 
to the conclusion. 



the explicit form of the energy density is specified in Eq. 
([3]) , the extension of the following formula to the general 
system is straightforward. 

Let us express Eq. in terms of the canonical 

variables because the following formula is based on the 
canonical theory developed in RefsJ^— . The magneti- 
zation dynamics without the spin torque and dissipation 
is described by the following Lagrangian densit y 28 ' 29 , 

M 

L= (cos 0-1) -E. (5) 

7 

The first term of Eq. (|S|) is the solid angle of the magne- 
tization dynamics in the spin space, or equivalently, the 
Berry phase. The canonical coordinate is q = tp, and the 
momentum is defined as p = dL /dq = — (M/7) (cos — 1). 
In terms of the canonical variables {q,p), the LLG equa- 
tion @ can be expressed as 

dq 1 dE aMH s d 

— — 11 1 ' II 

dt 1 + a 2 dp 1 + a 2 dp p 

ai _xgg 1 H S _ x d_ (6) 
(l + a 2 )M 9 dq l + a 29 dg m ' np ' 



II. FOKKER-PLANCK EQUATION FOR 
MAGNETIZATION 



25-27 



The LLG equation is given by=£r 
= —7m x H — jH s m x (rip x m) 



dm 

omx — , (2) 



where the gyromagnetic ratio and the Gilbert damping 
constant are denoted as 7 and a, respectively, m = 
(sin 9 cos tp, sin 9 sin tp, cos 9) and n p are the unit vectors 
pointing to the directions of the magnetization of the 
free and pinned layers, respectively. H = —dE/d(Mm) 
is the magnetic field, where M and E are the saturation 
magnetization and the magnetic energy density. For an 
in-plane magnetized system, E is given by 



E = ^ (m • e z f + 2irM 2 (m • e,) 2 

= _ MHl£ z 2 + 2i tM 2 (1 - z 2 ) cos 2 tp 



(3) 



where z — cos 9, and and 4ttM axe the uniaxial 
anisotropy field along the easy (z) axis and the demag- 
netization along the hard (x) axis, respectively. For the 
uniaxially anisotropic system discussed in Sec. IIV1 the 
demagnetization field should be absent. In Sees. IIVI and 
3 we assume that n p = e z . The strength of the spin 
torque in the unit of the magnetic field, H s , is given by 



H s = 



fvql 
2eMV ' 



(4) 



where r\ and V are the spin polarization of the current 
and the volume of the free layer, respectively. Although 



dp 
~dt 



1 dE aMH s d 



dq 1 



dq 



m • n R 



aM 



dE M 2 H S 



(l + a 2 )7' dp (1 + a 2 )j dp 



d 

g—mn 



(7) 



where g = sin 2 9. It should be noted that the spin torque 
term cannot be expressed as a gradient of the potential 
energy E, in general, and can be regarded as a damping 
or anti-damping. 

At a finite temperature, the random torque, —7m x 
h, should be added to the right-hand side of Eq. (J2J, 
where the components of the random field h satisfy the 
fluctuation-dissipation theorem^, 



{Whip)) 



2D 



S l3 S(t-t'). 



(8) 



Here D = a^k^T / (MV) is the diffusion coefficient. Due 
to the random torque, the magnetization switching can 
be regarded as the two-dimensional Brownian motion of 
a point particle in the (q,p) phase space. 

Let us introduce the distribution function of the mag- 
netization, W, in the (q,p) phase space. The Fokker- 
Planck equation is given by£ 



dW 
~dt 



dqdt 1 



dpdt 1 



D 



D 



or 



__) —g — W, 

7 / dp dp ' 



(9) 



where dq/ dt and dp/dt should be regarded as the right- 
hand sides of Eqs. and ([7]). In equilibrium (dW/dt = 
and H s = 0), the distribution function is identical to the 
Boltzmann distribution function (oc exp[— E/iksT)]). 
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III. WKB APPROXIMATION 



In general, the distribution function determined by Eq. 
([9|) depends on the two variables, (q,p), and it is very 
difficult to solve Eq. ^ for an arbitrary system. Thus, 
we use the following two assumptions. 

First, the low temperature assumption corresponding 
to the WKB approximations in RefsJ^r— is employed. 
We assume that the distribution function takes the fol- 
lowing for m 19 ' 21 , 



W oc exp {-aS/D) 



(10) 



In the zero temperature limit (a/D <C 1), Eq. Q is 
reduced to2£ 



dS_ 

dt 



dqdS_ _ _ x 
dt dq 



dpdS 



dS 



dt dp a9 \i ) \dp 



m\ 2 /as ,N 2 



(11) 



Here S and —dS/ dt = ffl can be regarded as the effective 
action and the Hamiltonian density, respectivel y 11 ' 21 " —. 
The corresponding Lagrangian density is then given by 



Jz? = — q\ q — p\ p — W, 



(12) 



dE 

' dq 



aMH s d a~i _j dE 

TTq^ 2 9p m ' np+ (1 + a 2 )M 9 ~dq~ 



1 + cr dq 



dE 
' dp 



aMH s d 

— — j — m-n p 

1 + cr dq 

M 2 H S d 
~(1 + a 2 )~j 9 dp~ 

2 



aM 



dE 



{l + a 2 )Y dp 



mn r 



+ A 



+ A 



-9 



dE 



M\ 2 fdE s 2 



dp 



(14) 



We perform the time average of j£f over 
the constant energy orbit, Jz? = J^dt^f/r = 

-[M/( 7 r)] J Q 27r dip^/ (dE/dz). Here we use the fact that 
the constant orbit is determined by the non-perturbative 
equation of motion, dip/dt — — (j/M)(dE/dz). The 
averaged Lagrangian density is expressed as 



M 



-EXi 



<7L } 



(15) 



where we renormalize Xe as {MXe M) — > Xe by which A e 
becomes a dimensionless variable. The effective Hamil- 
tonian is given by 



We = - F s + X E (1 + X E )^-^F a , 

1 + or 1 + or 

where F s and F a are given by 



2;r 



dp) 



dE/dz 



dE 



(1-^) ^ "« ^ 



dz 



dE 



d(p 



(16) 



(17) 



where X q = —dS/dq and A p = —dS/dp conjugated to q 
and p are the counting variables^ 2 -, 
is given by S = J dt% '. 



The effective action 



Second, we utilize the fact that the Gilbert damping 
constant and the spin torque strength are small. The 
values of the Gilbert damping constant a for the con- 
ventional ferromagnetic materials^! such as Co, Fe, and 
Ni are on the order of 10~ 2 . Also, the critical current 
of the spin torque switching 2 ^ is proportional to a|H|. 
Then, the switching time of the magnetization is much 
longer than the precession period t, and the energy E is 
approximately conserved during one period of the preces- 
sion. Following Refj 2 ^, we introduce the new canonical 
variables (E, s) accompanied by the new counting vari- 
ables (Xe, A s ). Then, the Lagrangian density is expressed 
as [see Appendix [A~] 



if 



dE 
~dt' 



1/0 



(13) 



-M 2 



T dip 2 )f — \ 1 (dE 

dE/dz ^ l ~ Z >\dz~) + l-z 2 \d(p 

The term proportional to H s in Eq. (|16p describes the 
work done by the spin torque, while the term propor- 
tional to cxXe in the second term of Eq. (|T^|) describes 
the energy dissipation due to the Gilbert damping [Ap- 
pendix [B]. The term proportional to aX 2 E arises from the 
thermal fluctuation. 

The switching barrier is obtained by integrating the 
Lagrangian density along the switching path. It is suffi- 
cient to calculate the integral along the optimal path 2 ^ in 
the low temperature limit, where the optimal path corre- 
sponds to the switching path obtained by Eq. ^ with- 
out the thermal fluctuation. The optimal path 2 ^ can be 
found by solving W e — 0. One of the solutions is given 
by E — — MHk/2, where z — ±1 and ip is arbitrary, 
due to which F s = F a = 0. The other two solutions of 
We — are given by Xe = 0, A^, where 



H S F S 
aMF a 



(19) 
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energy, 2E/(MHk) 



energy, 2E/(MHk) 



FIG. 1: (a) A typical dependence of \* E on the energy E for 
the uniaxially anisotropic system. The energy is normalized 
by the factor MHk/2- The dotted line corresponds to \ E = 0. 
The switching barrier is obtained by integrating \* E from -E m in 
to -Emax, i.e., the shaded region. The lower boundary -B m i n is 
fixed to —MHk/2, while the upper boundary -B ma x located at 
-(MH K /2)(I/I c f. (b) The dependence of \* E for the current 
I/I c = 0.2, 0.4, 0.6, 0.8 and 1.0. 

The equation of motion along Xe — describes the drift 
in the energy space due to the competition between the 
spin torque and the Gilbert damping at zero tempera- 
ture. The energy change E = —d-J^E/dXE along A^ = 
is given by 



E 



~f(H s F s -aMF a ) 



1 



(20) 



On the other hand, X* E corresponds to the time reversal 
path of Xe = 022, and thus, E\\ e= \* e = — E\x B =o- The 
switching barrier, A = aS/D, is then given by 

rB m « 



A 



V 



dEX% 



(21) 



The boundaries of the integral in Eq. (|21|) are deter- 
mined as follows. According to its definition, Eq. (fT^jl . 
A^ depends on the work done by spin torque and the en- 
ergy dissipation due to the damping, see also Appendix 
iBl In the region A^ > 0, the spin torque overcomes the 
damping, and the magnetization can move from the ini- 
tial state parallel to the easy axis without thermal fluc- 
tuation. On the other hand, the damping exceeds the 
spin torque in the region A^ < 0, and thus, the thermal 
fluctuation is required to achieve the switching. Then, 
the integral in Eq. (|2"Tj) in the region A^ < gives the 
switching barrier, i.e., the boundaries of the integral in 
Eq. (f21~T) are those of X* E < 0. In the latter sections, 
examples of X* E are shown [see Figs. [I]and|S]. 

Equation ([2Tj) with Eqs. (fT9)) and (|20l) is the main re- 
sult in this section, and enables us to calculate the current 
dependence of the switching barrier in the next sections. 
Equation (f2"Tj) is similar to Eq. (19) of Ref except for 
the additional condition on the integral range determined 
by Eq. ([2TJ|) . We emphasize that this difference is crucial 
to obtain the exponent b, as shown in Sees. IIVI and fVl 



IV. UNIAXIALLY ANISOTROPIC SYSTEM 

In this section, we calculate the switching barrier 
of the uniaxially anisotropic system, in which the en- 
ergy density is given by E = — (MH^/2) cos 2 9. Here 




FIG. 2: The dependence of the switching barrier of the 
uniaxially anisotropic system on the current I/I c . The barrier 
height is normalized by Ao. 

z = cos 8 can be expressed in terms of E and Hk as 
z = \] — IE I [MHk). The metastable states of the mag- 
netization locate at m = ±e 2 , and the initial state is 
taken to be m = e z . The functions F s and F a (Eqs. (fl7|) 
and (|18[) ) of this system are given by 



F s = 7 H K z (1 - 



F n = 



M 



z 2 (l-z 2 ), 



(22) 



(23) 



respectively. The critical current I c of the uniaxially 
anisotropic system is given by 



2aeMV 
ic — — r Hk- 



(24) 



The switching path, X* E in Eq. (flT))) , is given by 
1 + [I/{I c z)\. Figure Q] (a) shows a typical de- 
pendence of X* E on the energy E. The region A^ < is 
from E min = -MH K /2 to £ max = -(MH K /2)(I / I c f . 
Figure Q] (b) shows A^ for various currents. The upper 
boundary of A^ < 0, £ max , is zero at I = 0, and ap- 
proaches to —MHk/2 with increasing current. At I = I c , 
£ max equals -Emin. This behavior is particular for the 
uniaxially anisotropic system, where the switching can 
be reduced to the one-dimensional problem due to the 
rotational symmetry along the z axis. In this case, the 
effect of the spin torque can be described by the effective 
potential^ £ eff = -(MH K /2) cos 2 6 + (MH s /a)cos6. 
The effective potential has two minima at 8 — 0, tt and 
one maximum at 9 — cos _1 (/// c ). 

The switching barrier, Eq. (|2~T|) . is obtained by inte- 
grating the region of A^ < 0, i.e., the shaded region of 
Fig. Q] (a), and is given by 



A = A 1 



(25) 



where A = M Hk/ '(2A;bT) is the thermal stability. Fig- 
ure [2] shows the dependence of the switching barrier on 
the current I/I c . The barrier height is normalized by A . 
The exponent of the term 1 — I / I c is b = 2 in this system, 
which is consistent with Refs<^^. 
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(a) 



(b) 



(a) 




hard axis 

easy axis 




90 180 
angle from easy axis, 6 (deg) 
FIG. 3: (a) A typical energy map of an in-plane magnetized 
system. There are two minima (E — —MHk/2) at 9 — 0,n, 
two maxima (E = 2-kM 2 ) at (0, ip) = (tt/2, 0), (tt/2, 7r), and 
two saddle points (E = 0) at (0,cp) = (tt/2, tt/2), (tt/2, 3tt/2). 
The blue and red regions correspond to the stable (E < 0) 
and unstable (E > 0) regions, respectively, (b) A typical 
switching orbit in the in-plane magnetized system. The value 
of the parameters are M = 1000 emu/c.c, Hk = 200 Oe, 
a = 0.01, and T = 20 K, see Ref.— . 



V. IN-PLANE MAGNETIZED SYSTEM 

In this section, we investigate the switching barrier of 
an in-plane magnetized system. The energy density of 
this system is given by Eq. ([3]). The magnitudes of 
the demagnetization and the uniaxial anisotropy fields 
of the ferromagnetic materials of conventional in-plane 
Spin RAM are on the order of 1 T and 100 Oe£, re- 
spectively. Thus, Hyi/ (A-kM) is on the order of 10 ~ 2 . 
Figure [3] (a) shows a typical energy map of an in-plane 
magnetized system. There are two minima of the en- 
ergy (E = —MHk/2) at 9 = 0,tt and two maxima 
(E = 2nM 2 ) at (0,<p) = (tt/2, 0), (tt/2, tt). Because 
of the large demagnetization field, the magnetization 
switches through one of the saddle points (E = 0) at 
(0,<p) = (tt/2, tt/2), (tt/2, 3tt/2). Figure [3J (b) shows a 
typical switching orbit obtained by numerically solving 
the LLG equation^ 2 -. Starting from m = e z , the mag- 
netization precesses around the easy axis, and gradually 
approaches the saddle point. Then, the magnetization 
passes close to the saddle point, and relaxes to the other 
stable state, m = — e z . Since the switching time is dom- 
inated by the time spent during the precession around 
the easy axis in which —MHk/2 < E < 0, it is sufficient 
to evaluate \* E of Eq. (|2"Tj) in this energy region. 

The details of the calculation of the switching barrier 
are as follows. The variable z relates to E, Hk, and A-kM 



UnM cos 2 ip -2E/M 
AttM cos 2 (p + H K 



(26) 



Then, the functions F s and F a (Eqs. JXTJ) and ([Tgj)) arc 
given by 



2tt(H k + 2E/M) 
Ty/H K (H K +4nMy 



(27) 




0.9998, 

(b) 
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(d) 
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10 20 30 40 50 
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FIG. 4: The magnetization dynamics at zero temperature. 
The time evolution of m z (left) and the dynamic orbit (right) 
are shown, (a) At / = I c . The magnetization oscillates 
around the initial state (m = e z ) with small amplitude, (b), 

(c) At I = r(Ic - I )+Ic (r = 0.4 for (b) and r = 0.8 for (c)). 
The oscillation amplitude increases with increasing current. 

(d) At I = I*, the switching occurs. 



4 UirM-2E/M 
Wkf V 



AttM(H k + 2E/M)\ 
H k (4itM-2E/M) J 

4:nM(H K + 2E/M)\ 
H k (4ttM-2E/M) J 



where the precession period is given by 




(28) 



j^H K (4irM-2E/M) 



UirM(H K + 2E/M)\ 
H k (AttM-2E/M) J ' 
(29) 



The complete elliptic integrals of the first and second 
kinds are defined as K(fc) = J" dy/y/ (l — y 2 )(l — k 2 y 2 ) 

and E(fc) = J Q dy^f (1 — k 2 y 2 )/(l — y 2 ), respectively. In 
the limit of AttM — > 0, r, F s , and F a are identical to 
those calculated for the uniaxially anisotropic system. 



It should be noted that X* E of Eq. (|19p satisfies the 
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(a) 



(b) 






Emin 


Emax 


A E =0 






(lc<Klc*) 







energy, 2E/(MH K ) 



(d) 



r=0.2, 0.4, 0.6, 0.8, 1.0 (l c <l<l c ) 



=0.2, 0.4, 0.6, 0.8, 1.0{l<l c ) 



energy, 2E/(MH K ) 



energy, 2E/(MH K ) 



FIG. 5: (a) Typical dependence of X E of the in-plane mag- 
netized system on the energy E for I < I c . The switching 
barrier is obtained by integrating X* E from E m - m = —MHk/2 
to -Emax = 0, i.e., the shaded region, (b) Typical dependence 
of A * E on the energy for I c < I < I* . The lower boundary 
of the integral, E m i n , locates at —MHk/2 < 75 m in < 0. (c) 
X E for I = 7 C and 1 = 1*. (d) The dependence of \* E on the 
various currents, 7 = rl c for I < 7 C and I = r(7* — J c ) + 7 C 
for I c < I < I* where r = 0.2, 0.4, 0.6, 0.8 and 1.0. 



following relations: 



lim 

-MH K /2 



lim Ap = - i 1 - — 

£->0 



(30) 



(31) 



where the critical currents 7 C and I* are, respectively, 
given by 



2aeMV fTT „ , ^ 

I c = — r H K + 2ttM . 

nr] 



AaeMV 



y/AwM(H K +4ttM). 



(32) 



(33) 



Below, we concentrate on the region of 7 C < J* 
(Hit/ (in M) < 0.196). In the conventional ferromagnetic 
thin film, 77k and 47rM are on the order of 100 Oe and 
1 T, and thus, this condition is usually satisfied. For an 
infinite demagnetization field limit, I*/I c — 4/7r = 1.27... 

The physical meanings of 7 C and 7* are as follows. In 
Fig. 4, we show the magnetization dynamics at the zero 
temperature. The parameters are M = 1000 emu/c.c, 
H K = 200 Oe, 7 = 17.64 MHz/Oe, a = 0.01, r) = 0.8, 
and V = 7T x 80 x 35 x 2.5 nm 3 , respectively-^, by which 
7 C = 0.54 mA and I* — 0.68 mA, respectively. For I < I c , 
the damping overcomes the spin torque, and the initial 
state parallel to the easy axis is stable. At I — I c , the 
initial state becomes unstable, and the magnetization os- 
cillates around the easy axis with a small constant ampli- 
tude [see Fig. 2] (a)]. This has been already pointed out 




0.4 0.6 0.8 1.0 

current, l/l c * 

FIG. 6: The current dependence of the switching barrier 
normalized by Ao, i.e., A/An,. The current magnitude is nor- 
malized by 7*. 

in Reff2^, and I c has been considered to be the critical 
current of the magnetization switching. However, we em- 
phasize that I — I c is not a critical point, because I = I c 
does not guarantee the switching—. At I c < I < I*, the 
oscillation amplitude increases with increasing current, 
as shown in Figs. 0] (b) and (c). The critical current is 
7* , over which the magnetization can switch its direction 
without the thermal fluctuation, as shown in Fig. U](d). 

In Fig. [SJ we show the dependence of on the energy 
for various currents. A typical X* E for I < I c is shown 
in Fig. [S] (a). Here, X* E is always negative because the 
damping exceeds the spin torque. The switching barrier 
is obtained by integrating from 7J m i n = —MHk/2 
to £ max = 0. On the other hand, Fig. [S] (b) shows a 
typical \* E for 7 C < I < I*. From the initial state (E — 
—MHjt/2) to a certain energy Ti'min, the magnetization 
can move without the thermal fluctuation because the 
spin torque overcomes the damping (X E > 0). At E = 
E m i n , the magnetization dynamics is on a stable orbit, 
where the spin torque and the damping are balanced. 
The thermal fluctuation is required from E m i n to E max 
to switch the magnetization direction. In Fig. [S](c), we 
7 C and I = 7 C *. At I = I c , \* E < 0, 
-MH K /2. On the other hand, 
at I = I*, \* E > 0, and X* E = at E = 0. Figure[5](d) 
shows A^ for various currents. 

Figure [5] shows the current dependence of the switch- 
ing barrier obtained by numerically integrating the X E 
shown in Fig. [5] (d) [see Eq. (|21l) ]. in which the barrier 
height is normalized by Ao- The upper boundary of the 
integral range, E max , is taken to be E — 0. On the other 
hand, the lower boundary, E m i n , is —MH-k/2 for I < 7 C , 
while it is determined by numerically solving \ E = 
for I c < I < I*. We find that the current dependence 
of the switching barrier is well described by (1 — I/I*) b 
[see also Appendix [C] . The dependence is approximately 
linear (b ~ 1) for 7 < 7 C showing the consistence with 
Refsj 15 ' 16 . On the other hand, we find a nonlinear de- 
pendence for 7 C < 7 < 7* . The solid line in Fig. [7] shows 
the dependence of the exponent b on the current, where 
b is determined by the switching barrier shown in Fig. 
M as b = ln(A/A )/ln(l - 7/1*). As shown, b slightly 
increases with increasing current for 7 < 7 C (~ 0.87*). 
For 7 C < 7 < 7*, b rapidly increases, and reaches b > 2 



show X* E for 7 
and X* E is zero at E 



7 




0.06 



FIG. 7: The estimated value of the exponent b by analyzing 
Fig. ©with a function ln(A/A )/ ln(l - I/I*). The values of 
Hk in the solid, dashed, and dotted lines are 200, 500, and 
2250 Oe, respectively. 

near I < I*. It should also be noted that b is not uni- 
versal. The values of b with different Hk values are also 
shown in Fig. [7] by the dashed (Hk = 500 Oe) and dot- 
ted (H K = 2250 Oe ~ 0.18 x AirM) lines. As shown, 
the value of b decreases with increasing Hk- A further 
increase of Hk breaks the condition (J c < /*, or equiv- 
alently, Hk/(4itM) < 0.196), and beyond the scope of 
this paper. 

The nonlinear dependence of the switching barrier is 
important for an accurate estimation of the thermal sta- 
bility. Typical experiments are performed in the large 
current region, / < I*, to quickly measure the switch- 
ing. By applying a linear fit to Fig. [71 as done in the 
analysis of experiments^, the thermal stability would be 
significantly underestimated. 



VI. COMPARISON WITH PREVIOUS WORKS 

In Sec. IIII1 we assume that the switching time from the 
initial state (E = —MHk/2) to the saddle point (E = 0) 
is much longer than the precession period r. Since the 
precession period (Eq. ((29)) ) diverges at E = 0, this 
assumption is apparently violated in the vicinity of the 
saddle point. The condition under which our approxima- 
tion is valid is 



^l| E =o« 



MH K 



(34) 



By using Eqs. ®, (J27J, and 0, we find that Eq. ((MJ) 
can be expressed as 



AttM < 



H K 
16a 2 



(35) 



The parameters used in Sec. [V] satisfy this condition. 

According to Eq. (1351) . the present formula does 
not work for an infinite demagnetization field limit, 
AttM/Hk — > oo, where the switching occurs completely 
in-plane without precession. Then, the switching bar- 
rier is given by Aq(1 — I/I c ) 2 with 6 = 2 and I c = 
2aeMVHK/ (fir]), as shown in our previous work^. Thus, 
the present work is valid for Eq. (|35[) . while the previous 
work is valid for AttM/Hk —> oo. The switching barrier 



< 
5 



0.04 



0.02 




0.92 



0.94 0.96 

current, l/l c * 



0.98 0.99 



FIG. 8: The dots are the current dependence of the switching 
barrier in 0.92/,* < I < 0.99/*. The solid line represents a 
fitting proportional to (1 — I/Ic) b , where b' is assumed to be 
constant in this current region. 

in the intermediate region, Hk/ (16a 2 ) < AirM, remains 
unclear, and is beyond the scope of this paper. 

Let us also discuss the relation between the present 
work and one of our previous works in Ref^. In Ref»e2, 
we estimate that b ~ 3 for the in-plane magnetized sys- 
tem by numerically solving the LLG equation at various 
temperatures and adopting the phenomenological model 
of the switching^. The parameters are identical to those 
used in Fig. [Jj The switching currents at 0.1 K and 
20 K are 0.67 mA (~ 0.99/ c *) and 0.62 mA (~ 0.927*), 
respectively^. 

The main difference between the present work and 
Refs i 12 ' 17 is that, in Rets . 12 ' 17 the exponent b is assumed 
to be constant. In Fig. [8l the switching barrier shown in 
Fig. [6]is fitted by a function proportional to (1 — I / 'I*) b , 
with a constant b'. The current range is from 0.927* 
to 0.99/*, according to Refji 2 -. The obtained value of 
b' is 2.5, which is close to the estimated value (b ~ 3) 
in Refi^2. The difference may arise from the tempera- 
ture dependence of A, or the current dependence of the 
attempt frequency neglected in the present work. 



VII. SUMMARY 

In summary, we have developed a theory of the spin 
torque switching of the in-plane magnetized system based 
on the Fokker-Planck theory with WKB approximation. 
We derived the analytical expressions of the critical cur- 
rents, / c and I*. The initial state parallel to the easy 
axis becomes unstable at / = / c , which has been derived 
in RefM. On the other hand, at / = J*(~ 1.27/ c ), the 
switching occurs without the thermal fluctuation. Wc 
also find that the current dependence of the switching 
barrier is well described by (1 — I/I*) b , where the value 
of the exponent b is approximately linear for the cur- 
rent / < / c , while b rapidly increases with increasing 
current for / c < / < /*. The nonlinear dependence for 
/ c < / < I* is important for an accurate evaluation of 
the thermal stability because most experiments are per- 
formed in the current region of I c < / < I*. 
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Appendix A: Derivation of Eqs. (|13f) and (|14p 

First, we divide the Lagrangian density Jzf, Eq. ([T2"|) . 
into non-perturbative (-S?o) an d perturbative (Jz?i) parts. 
By using the explicit forms of dq/dt and dp/ dt in (Eqs. 
© and (J7J)), the non-perturbative Lagrangian density is 
given by 



dq 1 dE 
dt 1 + a 2 dp 



1 dE 



dp 

dt^l + a 2 ~dq 



(Al) 



On the other hand, the perturbative Lagrangian density 
is given by 



Jz?i — — A q 



+ A p 



aMH s d 
1 + a 2 dp 



m-rir, 



a-, 



_ x dE 



{l + a 2 )M- 



iH s _ x d 

1 + a 2 dq 
aMH, d 



aM 



1 + a 2 d 
M 2 H S 



dE 



(1 + a 2 )j dp 



(1 + a 2 )7 dp 



d 

5^-m-n p 



\2 " 

H 1 + a 2 



A 



a 



1 + a 2 



fJ 



M 



(A2) 



The canonical transformation from (q,p) to {E,s) 
is accompanied by the canonical transformation of the 
counting variables as X q = (dE/dq)\E + (ds/dq)X s and 
A p = (dE/dp)\E + {ds/dp)\ s . The non-perturbative 
Lagrangian is then given by J?fo = —\E(dE/dt) — 
X s [(ds/dt) - [s,H] P ], where d/dt = (d/dq)(dq/dt) + 
(d I 'dp) (dp I 'dt) and [, ]p is the Poisson bracket. In the 
small damping limit, s can be regarded as a physical time 
t. Since s is a conjugated variable of E ([s,H]p = 1), 
the non-perturbative Lagrangian density is given by Eq. 
(1131) . Since the switching barrier is determined by Jz?o 
and is independent of A s , we set A s = 0, according to 
ReL22. Then, J*f E = -%\ is given by Eq. fTgj). 



Appendix B: Work done by spin torque and 
damping 

Here we show that the functions F s and F a are propor- 
tional to the work done by the spin torque and damping. 



The time evolution of the magnetic energy, dE/dt 
■MH • (dm/dt), is calculated by using Eq. @ as 



dE 

— = W S + W a , 
dt 

where W s and W a are given by 



(Bl) 



W s = 



[n p • H— (m • n p ) (m ■ H) — an p ■ (m x H)] , 

(B2) 



1 + a 2 



H (m Hy 



(B3) 



Here, W s is the work done by spin torque while W a is 
the energy dissipation due to the damping. It should 
be noted that W a is always negative, while W s is ei- 
ther positive or negative depending on the current di- 
rection. Let us consider the in-plane magnetized sys- 
tem as an example, where the magnetic field is given by 
H = (— 4TrMm x , 0, H^rriz). Then, the time averages of 
W s and Wq over one period of the precession around the 
easy axis are given by 



w. - i^s 



:E a , 



(B4) 



(B5) 



1- 

respectively, where F s and F a are given by Eqs. (|27l) and 
(|2"51) . respectively. Thus, we can verify that F s and F a 
are proportional to the work done by the spin torque and 
the damping, respectively. 



Appendix C: Current dependence of the switching 
barrier near H 



Here, we derive an analytical formula of the switch- 
ing barrier near I*. For simplicity, we use the nor- 
malized variables k — H K /4irM, e = E/4irM 2 , and 
s = Hs/AttM. First, let us derive the approximated form 
of Xe (Eq. IT5)) ) around e ~ 0. The Taylor expansions 
of H s F s t and aMF a r are respectively given by 



H s F s t 2n(k + 2e)s 



y/k{l + k) ' 



(CI) 



aMF a r 
4ttM 



2(1 - k)s 
7^ 



1 - In 



(1 + k)e 
8k 



(C2) 

where e In e term appears and thus it is non-analytic at 
e = 0. Then, the approximated A^ is given by 



Ai 



Ie 
2I*k 



3 + k+ (1 - fc)ln 



(l + k)e 



8k 



(C3) 
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The energy -E m ; n = 4irM 2 Eo corresponding to the inter- 
section of X* E and Xe = can be obtained by solving the 
following self-consistent relation: 

F = { 1 _L\ 2 3KB 

V IS J 3 + k + (1 - k) ln[-(l + k)e /(8k)} ' 

(C4) 

Up to the first order of E m { a , the switching barrier is 
given by 



where E m \ n — inJVPeo is determined by Eq. (|C4|) . Equa- 
tion (|C5|) is only valid close to the switching current /* 
where E m [ n locates near E = 0. 
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